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Although human coronavirus OC43-OC43 (HCoV-OC43) is the coronavirus most commonly associated with 
human infections, little is known about its molecular epidemiology and evolution. We conducted a molecular 
epidemiology study to investigate different genotypes and potential recombination in HCoV-OC43. Twenty-nine 
HCoV-OC43 strains from nasopharyngeal aspirates, collected from 2004 to 2011, were subjected to RNA-dependent 
RNA polymerase (RdRp), spike, and nucleocapsid gene analysis. Phylogenetic analysis showed at least three distinct 
clusters of HCoV-OC43, although 10 unusual strains displayed incongruent phylogenetic positions between RdRp 
and spike genes. This suggested the presence of four HCoV-OC43 genotypes (A to D), with genotype D most likely 
arising from recombination. The complete genome sequencing of two genotype C and D strains and bootscan 
analysis showed recombination events between genotypes B and C in the generation of genotype D. Of the 29 strains, 
none belonged to the more ancient genotype A, 5 from 2004 belonged to genotype B, 15 from 2004 to 2006 belonged 
to genotype C, and 1 from 2004 and all 8 from 2008 to 2011 belonged to the recombinant genotype D. Molecular 
clock analysis using spike and nucleocapsid genes dated the most recent common ancestor of all genotypes to the 
1950s, genotype B and C to the 1980s, genotype B to the 1990s, and genotype C to the late 1990s to early 2000s, while 
the recombinant genotype D strains were detected as early as 2004. This represents the first study to describe 
natural recombination in HCoV-OC43 and the evolution of different genotypes over time, leading to the emergence 
of novel genotype D, which is associated with pneumonia in our elderly population. 


Coronaviruses cause infections in a wide variety of animals, 
resulting in respiratory, enteric, hepatic, and neurological dis¬ 
eases of various levels of severity. Based on genotypic and 
serological characterization, coronaviruses traditionally were 
classified into three distinct groups, groups 1, 2, and 3 (4). 
Recently, the Coronavirus Study Group of the International 
Committee for Taxonomy of Viruses has renamed the tradi¬ 
tional group 1, 2, and 3 coronaviruses as Alphacoronavirus, 
Betacoronavirus, and Gammacoronavirus, respectively (http: 
//talk.ictvonline .org/media/p/1230.aspx). 

The recent severe acute respiratory syndrome (SARS) epi¬ 
demic due to SARS coronavirus (SARS-CoV) and the identi¬ 
fication of SARS-related coronaviruses (SARSr-CoVs) from 
Himalayan palm civets and horseshoe bats in mainland China 
have led to a boost in interest in the study of coronaviruses in 
both humans and animals (5, 13, 24, 26, 33, 37, 55). Before the 
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SARS epidemic in 2003, there were only 19 known coronavi¬ 
ruses, including 2 human, 13 mammalian, and 4 avian corona¬ 
viruses. After the SARS epidemic, more than 20 additional 
novel coronaviruses have been described with complete ge¬ 
nome sequences (9, 24-26, 31, 42, 45, 50, 53, 54, 57). These 
include 3 human coronaviruses, 15 mammalian coronaviruses, 
and 4 avian coronaviruses. For human coronaviruses, human 
coronavirus NL63 (HCoV-NL63) (an alphacoronavirus) and 
human coronavirus HKU1 (HCoV-HKUl) (a betacoronavi¬ 
rus) have been discovered in addition to the two previously 
known human coronaviruses, human coronavirus 229E 
(HCoV-229E) (an alphacoronavirus) and human coronavirus 
OC43 (HCoV-OC43) (a betacoronavirus), as well as SARS- 
CoV (a betacoronavirus) (9, 45, 53, 56). While HCoV-229E 
and HCoV-OC43 were thought to account for 5 to 30% of 
human respiratory tract infections, HCoV-NL63 and HCoV- 
HKU1 often were detected in <5% of respiratory tract sam¬ 
ples (23, 29, 38). Outbreaks due to HCoV-OC43 also have 
been reported (3, 32, 44). Nevertheless, the different HCoVs 
often cocirculate, with one or two HCoVs being predominant 
depending on the geographical area and year (8, 11, 19, 23). 

Coronaviruses are unique in having a high frequency of 
homologous RNA recombination, which is a result of random 
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TABLE 1. Primers used for PCR and sequencing of the complete RdRp, S, and N genes of the 29 HCoV-OC43 strains 


Primer name 


Primer direction, sequence (5'-3') 


Gene target 


LPW 1114 
LPW 1118 
LPW 5188 
LPW 3709 
LPW 3064 
LPW 2579 
LPW 3064 
LPW 2579 
LPW 1223 
LPW 1127 
LPW 1162 
LPW 1166 
LPW 1261 
LPW 2094 
LPW 2095 
LPW 2098 
LPW 1839 
LPW 1178 
LPW 1177 
LPW 1183 
LPW 1275 
LPW 1189 
LPW 6547 
LPW 6548 
LPW 1192 
LPW 1045 
LPW 1195 
LPW 1198 


Forward, CYGTTTGTATATATTGCCGC 
Reverse, TGATTATCAAGTGTTAAAAC 
Forward, GCCCCTAGTTAGCGCTACTGAGTT 
Reverse, ACTTAGGATAATCCCAACCCAT 
Forward, CTGGGATGATATGTTACGCCG 
Reverse, GTGTGTTGTGAACARAAYTCRTG 
Forward, CTGGGATGATATGTTACGCCG 
Reverse, GTGTGTTGTGAACARAAYTCRTG 
Forward, TAAGTGCCTTTCAACAGGT 
Reverse, KGCCTTTTGCGTTTCTGC 
Forward, CCYRTTTGTRTGTATGATCC 
Reverse, YGCATAAAAAGTACCACC 
Forward, CTRCTATARYTATAGGTAGT 
Reverse, GCCCAAATTACCCAATTGTAGG 
Forward, TGATGCTGCTAAGATATATGG 
Reverse, ATTCCGARATAGCAATGCTGG 
Forward, ATCTTTTGTATGATTCTAATGG 
Reverse, GACACCAAGMCCATTAAT 
Forward, CWGCAGGTGTRCCATTTT 
Reverse, CCACAYTTCTTRAAACAAC 
Forward, TRAAATGGCCTTGGTATGT 
Reverse, TKWMWAGGAACTCTACAATA 
Forward, CTTCAAAGAACTATGGCATT 
Reverse, GACTGCAAATAGCCCAAATT 
Forward, AACCCMGAAACAAACAAC 
Reverse, GCAAGAATGGGGAACTGTGG 
Forward, GAGAGGCCCTAATCAGAA 
Reverse, TYAACTTCATTCATTTACTA 


RdRp 

RdRp 

RdRp 

RdRp 

RdRp 

RdRp 

RdRp 

RdRp 

RdRp 

RdRp 

S 

S 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

N 

N 

N 

N 


template switching during RNA replication that is thought to 
be mediated by a copy-choice mechanism (28, 46). Their ten¬ 
dency for recombination and high mutation rates may allow 
them to adapt to new hosts and ecological niches. During our 
previous investigations on the molecular epidemiology of 
HCoV-HKUl, we documented the first evidence for natural 
recombination in coronavirus associated with human infection, 
resulting in the generation of different HCoV-HKUl geno¬ 
types (23, 52, 56). Since some strains of HCoV-HKUl were 
found to display incongruent phylogenetic relationships upon 
the analysis of the RNA-dependent RNA polymerase (RdRp), 
spike (S), and nucleocapsid (N) genes, recombination events 
were suspected and later confirmed with the complete genome 
sequencing of 22 strains of HCoV-HKUl and recombination 
analysis (52). Although HCoV-OC43 is thought to be the most 
commonly encountered human coronavirus, no similar molec¬ 
ular epidemiology studies have been performed, and little is 
known about its evolution among humans. Only five complete 
genome sequences of HCoV-OC43, two from the same Amer¬ 
ican Type Culture Collection (ATCC) strain, VR759, that was 
isolated in 1967, one Paris strain that was isolated in 2001, and 
two Belgium strains detected in 2003 and 2004, were available 
in GenBank (39, 47, 48). In this study, we investigate the 
presence of different genotypes among HCoV-OC43 strains 
and identify potential recombination events that lead to the 
generation of novel genotypes, a situation analogous to that 
observed for HCoV-HKUl. HCoV-OC43 detected from the 
nasopharyngeal aspirates (NPAs) from 29 patients with respi¬ 
ratory tract infections from 2004 to 2011 were subjected to 
complete RdRp, S, and N gene sequencing and analysis. The 
clinical characteristics of patients also were analyzed in rela¬ 


tion to molecular epidemiology results. As initial analyses 
showed the presence of potential recombination events, two 
complete genomes of HCoV-OC43 were selected for sequenc¬ 
ing and further analysis. The emergence of a novel genotype of 
HCoV-OC43 through recombination and the evolution of dif¬ 
ferent HCoV-OC43 genotypes also was described. 

MATERIALS AND METHODS 

HCoV-OC43 strains. Twenty-nine HCoV-OC43 strains were detected from 
nasopharyngeal aspirates of patients with acute respiratory tract infections who 
had been admitted to two Hong Kong public hospitals, Queen Mary Hospital and 
Pamela Youde Nethersole Eastern Hospital. They were randomly selected dur¬ 
ing a 7-year period (November 2004 to February 2011) and were included in this 
study (23). All NPAs were negative for influenza A and B viruses, parainfluenza 
virus types 1, 2, and 3, respiratory syncytial virus, and adenovirus by direct 
immunofluorescence, and they also were negative for metapneumovirus, human 
coronavirus HKU1, human coronavirus 229E, and human coronavirus NL63 by 
reverse transcription-PCR (RT-PCR) (23, 56). 

RNA extraction. Viral RNA was extracted from the NPAs of the correspond¬ 
ing patients using a QIAamp viral RNA Minikit (QIAgen, Hilden, Germany). 
The RNA pellet was eluted in 50 |xl of DNase-free, RNase-free double-distilled 
water and was used as the template for RT-PCR. 

RT-PCR and sequencing of the complete RdRp, S, and N genes of HCoV- 
OC43 and phylogenetic analysis. The RNA was converted to cDNA by a com¬ 
bined random priming and oligo(dT) priming strategy. The complete RdRp, S, 
and N genes of HCoV-OC43 from 29 NPAs were amplified and sequenced using 
the primers shown in Table 1 and the strategy described in our previous publi¬ 
cations (23, 53). The nucleotide and the deduced amino acid sequences of the 
RdRp, S, and N genes were compared to those of HCoV-OC43 and other group 
2 coronaviruses. A phylogenetic tree was constructed using the maximum-like¬ 
lihood method in PhyML with GTR for RdRp and N genes and GTR+I for the 
S gene. The substitution models were selected according to Akaike information 
criterion (AIC) implemented in ModelGenerator version 0.85 (17). The robust¬ 
ness of branches was assessed by bootstrap analysis with 1,000 replicates. 
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TABLE 2. Calculated tMRCA values for S and N genes 


Genotype and 
analysis model 


Strict clock 


Relaxed uncorrelated clock 0 




Log-normal 


Exponential 

Mean 

HPD 

Mean 

HPD 

Mean 

HPD 

Constant size 

S gene 







A, B, and C 

48.1555 

44.7545-53.5444 

49.414 

44.4431-58.0039 

54.284 

44.5851-70.9238 

B and C 

27.4725 

22.2275-33.6636 

21.3901 

10.6186-33.4467 

26.9695 

14.7718^12.409 

A 

45.5514 

44.6245-46.6342 

45.62 

44.2208-47 

46.529 

44.3516-50.3417 

B 

11.4415 

9.5698-13.5421 

10.7938 

8.2477-14.2259 

14.8138 

10.0011-21.2453 

C 

8.121 

7.2977-9.0669 

8.2939 

7.1553-9.9183 

9.2664 

7.2557-12.4169 

N gene 







A, B and C 

55.6007 

44.102-80.1591 

55.2781 

44.099-79.5748 

54.433 

44.0908-77.5797 

B and C 

26.2353 

13.0594-44.625 

24.6398 

10.2003^13.6165 

21.6914 

9.1696^10.8566 

A 

45.58 

44.08-48.971 

45.5361 

44.08-48.9703 

45.5039 

44.08^19.1539 

B 

16.5939 

8.8253-27.7325 

15.3999 

8.2083-26.293 

13.2508 

8.1759-22.8651 

C 

13.6795 

8.2804-21.6243 

13.0869 

7.6379-21.0773 

12.1722 

7.4932-20.213 

Bayesian skyline 

S gene 







A, B and C 

45.0589 

44.4048-45.8426 

44.982 

44.09-46.3717 

45.0131 

44.1592^16.4873 

B and C 

20.8376 

17.939-23.7088 

10.1779 

8.5117-12.5169 

10.5825 

8.8659-12.8411 

A 

44.8137 

44.3591-45.3259 

44.4406 

44.1038-44.9446 

44.4342 

44.1312^14.8766 

B 

9.5314 

8.6617-10.4592 

8.4944 

8.0953-9.0724 

8.479 

8.1131-9.005 

C 

7.3604 

7.0823-7.6852 

7.3576 

7.08-7.8177 

7.3233 

7.08-7.7445 

N gene 

A, B, and C 

52.9938 

44.081-78.656 

52.6201 

44.0807-77.0595 

51.413 

44.0802-72.6384 

B and C 

23.9449 

12.055-42.3612 

23.1098 

10.0302-41.162 

20.571 

8.7883-38.149 

A 

45.151 

44.08-47.9311 

45.1336 

44.08-47.8964 

45.119 

44.08-48.0281 

B 

15.6571 

8.8518-25.9166 

15.1134 

8.3088-24.9999 

13.5672 

8.1754-23.0218 

C 

13.6085 

8.2093-21.2408 

13.3621 

7.8461-20.9611 

12.7738 

7.522-21.2622 


a The adopted model is shaded gray. 


Complete genome sequencing. Two complete genomes of HCoV-OC43 from 
two patients (HK04-01 and HK04-02) were amplified and sequenced using the RNA 
extracted from the original NPAs as templates. The RNA was converted to cDNA 
by a combined random priming and oligo(dT) priming strategy. The cDNA was 
amplified by degenerate primers designed by multiple alignments of available 
HCoV-OC43 genome sequences using strategies described in our previous publica¬ 
tions (24, 50, 52, 53), and the coronavirus database (CoVDB) (15) was used for 
sequence retrieval. These primer sequences are available on request. The 5' ends of 
the viral genomes were confirmed by the rapid amplification of cDNA ends (RACE) 
using a 573' RACE kit (Roche, Germany). Sequences were assembled and manually 
edited to produce final sequences of the viral genomes. 

Genome analysis. The nucleotide sequences of the genomes and the deduced 
amino acid sequences of the open reading frames (ORFs) were compared to 
those of other coronaviruses using CoVDB (15). Phylogenetic tree construction 
was performed using the neighbor-joining method using ClustalX 1.83, with 
bootstrap values calculated from 1,000 trees. The prediction of the receptor 
binding domain of spike protein was performed using InterProScan (2). The 
prediction of potential N-glycosylation sites in the spike proteins was performed 
using the CBS NetNGlyc 1.0 server (http://www.cbs.dtu.dk/services/NetNGlyc/). 
The number of synonymous substitutions per synonymous site, K s , and the 
number of nonsynonymous substitutions per nonsynonymous site, K a , for each 
coding region was calculated using the Nei and Gojobori substitution model with 
Jukes-Cantor correction in MEGA 4.0 (41). Bootscan analysis was used to detect 
possible recombination using the nucleotide alignment of the genome sequences 
of HCoV-OC43 strains HK04-01 and HK04-02, the ATCC strain, and Belgium 
strain BE03, generated by ClustalX version 1.83, and edited manually. Bootscan 
analysis was performed using Simplot version 3.5.1 as described previously (21, 
22, 52), with strains HK04-02 and BE04 as the query. 

Estimation of divergence dates. Divergence times for the three OC43 geno¬ 
types, A, B, and C, were calculated using a Bayesian Markov chain Monte Carlo 
(MCMC) approach as implemented in BEAST (version 1.6.1) as described 
previously (7, 22, 47, 49). One parametric model (constant size) and one non- 


parametric model (Bayesian skyline) tree priors were used for the inference. 
Analyses were performed under the GTR+I and GTR substitution models for S 
and N gene sequence data, respectively, and using both a strict and a relaxed 
molecular clock. For the S gene, the MCMC run was 3 X 10 7 steps long, while 
for the N gene the MCMC runs were 4 X 10 8 and 6 X 10 8 when the Bayesian 
skyline and constant-size coalescent tree prior was used, respectively, with sam¬ 
pling every 1,000 steps. Convergence was assessed on the basis of the effective 
sampling size after a 10% burn-in using Tracer software, version 1.5 (7). The 
mean time of the most recent common ancestor (tMRCA) and the highest 
posterior density regions at 95% (HPDs) were calculated, and the best-fitting 
models were selected by a Bayes factor using marginal likelihoods implemented 
in Tracer (Tables 2 and 3) (40). Constant population size under a relaxed-clock 
model with an uncorrelated exponential distribution was adopted for making 
inferences, as Bayes factor analysis for the S gene indicated that this model fitted 
the data better than other models tested (Table 3). The trees were summarized 
in a target tree by the Tree Annotator program included in the BEAST package 
by choosing the tree with the maximum sum of posterior probabilities (maximum 
clade credibility) after a 10% burn-in. 

Nucleotide sequence accession numbers. The nucleotide sequences of the two 
genomes of HCoV-OC43 have been lodged within the GenBank sequence da¬ 
tabase under accession no. JN129834 and JN129835. 


RESULTS 

Clinical characteristics of patients with HCoV-OC43 infec¬ 
tions. The clinical characteristics of the 29 patients with 
HCoV-OC43 infections are summarized in Table 4. Sixteen 
were males and 13 were females. Thirteen were children and 
16 were adults, with most patients at the extremes of age 
(median age, 36 years old; range, 24 days to 94 years old). 
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TABLE 3. Model selection results for the S and N genes of HCoV-OC43 through comparison of marginal likelihoods and log 10 Bayes factors 

for each pair of models 


Model 

combination 0 

Log-normal marginal likelihood 



Model combination 



P* (model/data) 

SE C 

Strict_CST 

StrictBSP 

Uced_CST 

UcedBSP 

Ucld_CST 

UcldBSP 

S gene 









Strict CST 

-7,420.992 

±0.119 


9.524 

-9.283 

-3.477 

-5.761 

-4.179 

Strict BSP 

-7,442.921 

±0.147 

-9.524 


-18.806 

-13 

-15.285 

-13.702 

Uced CST 

-7,399.618 

±0.143 

9.283 

18.806 

- 

5.806 

3.521 

5.104 

Uced BSP 

-7,412.987 

±0.155 

3.477 

13 

-5.806 


-2.285 

-0.702 

Ucld CST 

-7,407.726 

±0.169 

5.761 

15.285 

-3.521 

2.285 


1.583 

UcldBSP 

-7,411.371 

±0.167 

4.179 

13.702 

-5.104 

0.702 

-1.583 


N gene 









Strict CST 

-2,183.957 

±0.023 


-0.929 

-1.969 

-2.647 

-0.48 

-1.289 

Strict BSP 

-2,181.818 

±0.031 

0.929 


-1.04 

-1.718 

0.449 

-0.36 

Uced CST 

-2,179.423 

±0.03 

1.969 

1.04 


-0.678 

1.489 

0.681 

Uced BSP 

-2,177.861 

±0.028 

2.647 

1.718 

0.678 


2.167 

1.359 

Ucld CST 

-2,182.852 

±0.026 

0.48 

-0.449 

-1.489 

-2.167 


-0.809 

Ucld BSP 

-2,180.99 

±0.032 

1.289 

0.36 

-0.681 

-1.359 

0.809 


a Three molecular clock models, a strict clock and two relaxed clocks assuming either uncorrelated exponential (Uced) or uncorrelated lognormal distribution (Ucld) of 
substitution rates, were compared in combination with models of demographic history, constant size (CST), and Bayesian skyline (BSP). The adopted model is shaded gray. 

b Marginal likelihood estimated using the program Tracer 1.5. A log 10 Bayes factor was calculated for each pair of model combinations (i.e., model 1 in row versus 
model 2 in column). 

c Standard errors for the marginal likelihoods. 


Fifteen patients, especially infants and young children, pre- 36-year-old male and a 46-year-old female, respectively. The 

sented with upper respiratory tract infections. Of the 14 pa- 36-year-old male patient (patient 5) did not have underlying 

tients with pneumonia, three were infants or young children, disease and presented with chest pain and hemoptysis compli- 

nine were elderly (>60 years old), and the other two were a cated by right pleural effusion. The 46-year-old female patient 

TABLE 4. Clinical characteristics of the 29 patients with HCoV-OC43 infections 

Patient no. 

Mo/yr 

Strain no. 

Sex° 

Age 

Diagnosis* 

Genotype 

i 

11/2004 

HK04-01 

F 

9 yr 

URTI 

c 

2 

11/2004 

HK04-02 

M 

35 mo 

URTI 

D 

3 

11/2004 

HK04-03 

F 

1 yr 

URTI, febrile convulsion 

C 

4 

11/2004 

HK04-04 

M 

24 m 

URTI, febrile convulsion 

B 

5 

11/2004 

HK04-05 

M 

36 yr 

Pneumonia, pleural effusion 

C 

6 

11/2004 

HK04-06 

F 

82 yr 

Pneumonia, COPD exacerbation 

C 

7 

11/2004 

HK04-07 

F 

24 days 

URTI 

B 

8 

11/2004 

HK04-08 

F 

5 yr 

Pneumonia 

B 

9 

11/2004 

HK04-09 

M 

lyr 

URTI, febrile convulsion 

C 

10 

11/2004 

HK04-10 

M 

2 yr 

URTI 

B 

11 

11/2004 

HK04-11 

M 

85 yr 

URTI, urinary tract infection 

C 

12 

11/2004 

HK04-12 

M 

1 mo 

URTI 

C 

13 

11/2004 

HK04-13 

M 

76 yr 

URTI 

C 

14 

12/2004 

HK04-14 

F 

32 yr 

URTI 

C 

15 

12/2004 

HK04-15 

M 

72 yr 

URTI 

C 

16 

12/2004 

HK04-16 

M 

lyr 

Pneumonia 

C 

17 

12/2004 

HK04-17 

F 

46 yr 

Pneumonia, oral herpes 

B 

18 

12/2004 

HK04-18 

M 

68 yr 

Pneumonia, COPD exacerbation 

C 

19 

1/2005 

HK05-01 

M 

11 mo 

URTI, NSAID-induced angioedema 

C 

20 

12/2005 

HK05-02 

M 

88 yr 

Pneumonia 

C 

21 

1/2006 

HK06-01 

F 

84 yr 

URTI 

C 

22 

12/2008 

HK08-01 

F 

27 mo 

URTI 

D 

23 

11/2008 

HK08-02 

M 

83 yr 

Pneumonia, pleural effusion, CHF 

D 

24 

10/2009 

HK09-01 

M 

26 mo 

Pneumonia 

D 

25 

11/2009 

HK09-02 

M 

84 yr 

Pneumonia 

D 

26 

10/2010 

HK10-01 

F 

82 yr 

Pneumonia 

D 

27 

10/2010 

HK10-02 

F 

87 yr 

Pneumonia, acute pulmonary edema 

D 

28 

2/2011 

HK11-01 

F 

72 yr 

Pneumonia, CHF 

D 

29 

2/2011 

HK11-02 

F 

94 yr 

Pneumonia 

D 


a F, female; M, male. 

b URTI, upper respiratory tract infection; COPD, chronic obstructive pulmonary disease; NSA1D, nonsteroidal antiinflammatory drug; CHF, congestive heart failure. 
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RdRp 



. D 



N 


648 | HCoV-QC43-HK09-01 
I HCoV-OC43-HK11-02 

-HCoV-OC43-HK11-01 

HCoV-OC43-H K10-02 
HCoV-OC43-HK10-01 
HCoV-OC43-H K09-02 
HCoV-OC43-HK08-02 
HCoV-OC43-HK08-01 
HCoV-OC43-HK06-01 

-HCOV-OC43-HK05-02 

-HCoV-OC43-HK05-01 

-HCoV-OC43-HK04-18 

-HCoV-OC43-HK04-14 

640 1 -HCoV-OC43-BE04 (AY903460) 


' HCoV-OC43-HK04-02 

-HCoV-OC43-HK04-13 

-HCoV-OC43-HK04-09 

HCoV-OC43-HK04-11 
HCoV-OC43-HK04-06 
HCoV-OC43-HK04-03 
HCoV-OC43-HK04-01 
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FIG. 1. Phylogenetic analysis of the complete RdRp, S, and N genes of the 21 HCoV-OC43 strains. The trees were constructed by the 
maximum-likelihood method, and bootstrap values were calculated from 1,000 trees. A total of 2,783, 4,089, and 1,347 nucleotide positions in 
RdRp, S, and N, respectively, were included in the analysis. The scale bar indicates the estimated number of substitutions per 2,000, 1,000, or 20 
nucleotides as indicated. A, genotype/clade A; B, genotype/clade B; C, genotype/clade C; D, genotype D. Genotype D strains displaying 
incongruent phylogenetic relationships are in boldface. 


(patient 17) had a history of asthma, and her pneumonia was 
complicated by oral herpes. Two elderly patients were compli¬ 
cated by the exacerbation of chronic obstructive pulmonary 
diseases, and three were complicated by congestive heart fail¬ 
ure or acute pulmonary edema. Two children were compli¬ 
cated by febrile convulsion. Two patients (patients 5 and 18) 
had recent travel history to mainland China before symptom 
onset. Except for a 94-year-old female patient (patient 29) who 
died of pneumonia with superimposed Pseudomonas aerugi¬ 
nosa infection, no other respiratory pathogens were detected 
in other patients and all survived. 

RT-PCR and sequencing of the complete RNA-dependent 
RNA polymerase, spike, and nucleocapsid genes of HCoV- 
OC43 and phylogenetic analysis. The complete RdRp, S, and 
N genes of HCoV-OC43 from the 29 NPAs were amplified and 
sequenced. The analysis of the RdRp genes showed that the 29 
HK strains possessed 99.6 and 100% nucleotide identity to 
those of Belgium strains BE04 and BE03, respectively. Multi¬ 
ple alignments revealed that there were 33 nucleotide substi¬ 
tutions among the RdRp genes of the 29 HK strains. Phyloge¬ 
netic analysis showed the existence of three clusters of 
sequences, all with high bootstrap values of >900 (Fig. 1). One 
cluster, clade A, was formed by the ATCC and Paris strains. 
Fourteen HK strains were clustered with the BE03 and BE04 
strains, forming clade B. The other 19 HK strains were clus¬ 
tered together, forming clade C. 

The analysis of the S genes showed that 5 of the 29 HK 
strains from 2004 (HK04-04, HK04-07, HK04-08, HK04-10, 


and HK04-17) possessed high nucleotide identities (99.1 to 
99.5%) to Belgium strain BE03, while the other 24 sequences 
possessed high nucleotide identities (99.2 to 99.9%) to Bel¬ 
gium strain BE04. Multiple alignments revealed that there 
were 166 nucleotide substitutions among the S genes of the 29 
HK strains. Phylogenetic analysis showed that the sequences 
fell into three clusters, all with high bootstrap values of >990 
(Fig. 1). Again, one cluster, clade A, was formed by the ATCC 
and Paris strains. In line with results from pairwise identities 
described above, five HK strains from 2004 and Belgium strain 
BE03 formed another cluster, clade B, while the other 24 HK 
strains and Belgium strain BE04 formed the third cluster, 
clade C. 

Analysis of the N genes of the 29 HK strains showed results 
similar to those for the S genes. Multiple alignments revealed 
that there were 23 nucleotide substitutions among the N genes 
of the 29 HK strains. Phylogenetic analysis showed the exis¬ 
tence of three clusters of sequences (Fig. 1). Again, the ATCC 
and Paris strains formed a distinct cluster, clade A. Similarly to 
findings of S gene analysis, five HK strains from 2004 were 
most closely related to the Belgium strain BE03, forming clade 
B. The N gene sequences of the other 24 HK strains were more 
closely related to those of Belgium strain BE04, forming an¬ 
other cluster, clade C. 

From these results, 10 unusual strains, including Belgium 
strain BE04 and nine HK strains (one from 2004 and eight 
from 2008 to 2011), were found to display incongruent phylo¬ 
genetic positions upon the analysis of their RdRp, S, and N 
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gene sequences. They belonged to clade B in the RdRp tree 
but belonged to clade C in the S and N trees. Multiple align¬ 
ment also revealed that these 10 unusual strains displayed 
higher nucleotide similarity to other clade B strains in their 
RdRp genes, but they displayed higher nucleotide similarity to 
clade C strains in their S and N genes. These results suggested 
the presence of at least three genotypes of HCoV-OC43: ge¬ 
notype A (comprising only the ATCC and Paris strains), ge¬ 
notype B (including Belgium strain BE03 and five HK strains 
from 2004), and genotype C (including 15 HK strains from 
2004 to 2006). Moreover, 10 unusual strains, including Belgium 
strain BE04 and nine HK strains, may represent an additional 
genotype, genotype D, which may have arisen from recombi¬ 
nation between genotypes B and C. Interestingly, while one 
(HK04-02) of the nine genotype D strains from Hong Kong 
was detected in 2004, all eight HK strains from 2008 to 2011 
included in this study belonged to genotype D, suggesting that 
this genotype is the predominant genotype in recent years in 
our population. Moreover, seven of these eight recent strains 
were associated with pneumonia, especially in the elderly, 
three of which were complicated by heart failure (Table 4). 

Complete genome analysis of two HCoV-OC43 strains, 
HK04-01 and HK04-02. To determine if genotype D strains 
have arisen from recombination between genotypes B and C, 
the complete genome sequences of one selected strain from 
genotype D, HK04-02, and one selected strain from genotype 
C, HK04-01, were determined (complete genome sequence 
was available for genotype B strain BE03). The sizes of the 
genomes of the two HCoV-OC43 strains, HKU04-01 and 
HKU04-02, ranged from 30,710 to 30,722 nucleotides, with a 
G+C content of 37%. Their genome organizations were typi¬ 
cal of Betacoronavinis (Fig. 2). 

Phylogenetic trees using the nucleotide sequences of puta¬ 
tive proteins and polypeptides (nspl to nspl6, NS2a, hemag¬ 
glutinin-esterase [HE], S, NS5a, envelope [E], membrane [M], 
and N) of the two present strains and five other HCoV-OC43 
complete genome sequences available (from ATCC, Paris, and 
Belgium strains) were constructed and are shown in Fig. 3. In 
most of the trees, the ATCC and Paris strains formed a distinct 
cluster separate from the two Belgium and two HK strains. In 
most nsp proteins except nsp5 and nspl5, strains HK04-02 and 
BE04 were more closely clustered with strain BE03 than with 
strain HK04-01 (in nspl to nsp4, nsp6, nsp8, and nspl2) or 
clustered closely with both strains (in nsp7, nsp9 to nspll, 
nspl3, nspl4, and nspl6). However, in HE, S, E, M, and N 
genes, strains HK04-02 and BE04 were more closely clustered 
with strain HK04-01 than with strain BE03. This suggested that 
the two strains have arisen from recombination. 

Results from bootscan analysis were in line with the obser¬ 
vations described above. From the 5' end of the genome to 
position 22500, bootscan analysis showed a number of possible 
recombination sites when the genome of strain HK04-02 was 
used as the query (Fig. 2). Upstream of position 15500, most of 
the region exhibits higher bootstrap support for the clustering 
of strain HK04-02 with BE03, except between positions 2500 
and 5000, where higher bootstrap support for clustering with 
strain HK04-01 was observed. From positions 17000 to 22500, 
most of the region exhibits higher bootstrap support for clus¬ 
tering between strains HK04-02 and HK04-01. From position 
22500 to the 3' end of the genome, no further recombination 


sites were detected, with high bootstrap support for clustering 
between strains HK04-02 and HK04-01 observed. These indi¬ 
cated that recombination has taken place between nucleotide 
positions 2500 and 5000 (corresponding to the nsp2/nsp3 junc¬ 
tion), between nucleotide positions 15500 and 17000 (corre¬ 
sponding to the nspl2/nspl3 junction), and/or around nucleo¬ 
tide position 22600 (corresponding to the NS2a/HE junction). 
Similar results were obtained when strain BE04 was subjected 
to bootscan analysis (Fig. 2). Among these three potential 
recombination sites, the one at the nspl2/nspl3 junction most 
likely has resulted in the incongruent phylogenetic positions 
observed in the RdRp, S, and N genes among the genotype D 
strains, as this is the junction that bridges between reciprocal 
clustering to genotype B and C strains observed in RdRp and 
S genes, respectively. 

Further analysis of the present two HCoV-OC43 genomes 
revealed the presence of a putative transcription regulatory 
sequence (TRS) motif, 5'-UC(U/C)AAAC-3', at the 3' end of 
the leader sequence, and it preceded each translated ORF 
except for those of the NS5a and E genes, as described for 
HCoV-OC43 genomes previously (39). For NS5a, it has been 
suggested that 5'-UCIJ U A AG-3' was the putative TRS in the 
ATCC strain (39). However, this sequence was absent from the 
two Belgium strains, BE03 and BE04, and the present two 
strains, HK04-01 and HK04-02, as a result of a 12-nucleotide 
deletion from the region (Fig. 4). Upstream of this deletion, a 
potential alternative TRS, 5' - U C U A G C A - 3', was identified 20 
bp upstream to the start codon. It has been found previously 
that sequences in a similar region, upstream of E of HCoV- 
NL63 and NS5a of HCoV-HKUl, were homologous to a frag¬ 
ment of the corresponding leader sequences, which may serve 
as a compensation mechanism for the absence of an optimal 
TRS (34, 35). In HCoV-OC43, a 39-bp sequence upstream of 
NS5a also was found to be homologous to its leader sequence 
(Fig. 4). Further experiments are required to determine if 
NS5a utilizes a TRS or another mechanism for translation. As 
in sialodacryoadenitis virus (SDAV) and mouse hepatitis virus 
(MHV), the putative E gene of HCoV-OC43 may share the 
same TRS with NS5a, suggesting that the translation of the E 
protein is cap independent, possibly via an internal ribosomal 
entry site (IRES) (16). In MHV, an IRES element, UUUUA 
UUCUUUUU, has been identified upstream of the initiation 
codon of the E protein (16). A stretch of 13 nucleotides, 
CUUUAUUACCUGG, is present at this position in HCoV- 
OC43 according to multiple alignment. Further experiments 
would determine if this or another nearby sequence acts as an 
IRES for the E protein of HCoV-OC43. 

The analysis of the predicted S proteins of the present 29 
HCoV-OC43 HK strains revealed a potential N-terminal sig¬ 
nal peptide of 14 to 17 amino acids by SignalP-HMM and 
SignalP-NN, respectively. A potential cleavage site located af¬ 
ter RRSRR that is identical to that of bovine CoV (BCoV) (1) 
and the two Belgium HCoV-OC43 strains BE03 and BE04 
(48), between residues 766 and 767, where S will be cleaved 
into SI and S2, was identified in the S proteins of all 29 strains. 
The residues RRSRG were observed upstream of the potential 
cleavage site in the S proteins of the ATCC and Paris strains. 
It has been suggested that the G-to-R amino acid change in the 
last position, leading to an RRSRR motif, leads to an in¬ 
creased cleavability compared to that of the ATCC prototype 
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-HCoV-OC43-ATCC-VR759-AY391 111 (Genotype A) 

_ HCoV-OC43-BE03-AY903459 (Genotype B) 

. HCoV-OC43-HK04-01 (Genotype C) 



HK04-02 



BE04 



FIG. 2. Genome organization and bootscan analysis of the HCoV-OC43 genomes. Bootscanning was conducted with Simplot version 3.5.1 (F84 
model; window size, 1,000 bp; step, 200 bp) on a gapless nucleotide alignment, which was generated with ClustalX with the genome sequences of 
strains HK04-02 (upper) and BE04 (lower) as the query sequences. 
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FIG. 4. Alignment of the leader sequence and the homologous sequence upstream of the NS5a ORF. Sequence homology between the 
sequences is marked by asterisks. The putative TRS is highlighted in gray. The TRS suggested by St.-Jean et al. (39) are in boldface, and the 
corresponding deletions in BE03, BE04, HK04-01, and HK04-02 are indicated by the box. The initiation codons of NS5a are underlined. 


HCoV-OC43 spike protein (48). The analysis of the S proteins 
of HCoV-OC43 also reveals different numbers of potential 
N-linked glycosylation sites in different strains, with six NXS 
and eight NXT sites in the ATCC and Paris strains, six NXS 
and nine NXT sites in strain BE03, and seven NXS and nine 
NXT sites in strains BE04, HK04-01, and HK04-02. Sialic acid 
was known to be the receptor for S protein binding in HCoV- 
OC43, although the receptor-binding domain is not well de¬ 
fined (18). The S gene of the 29 HK strains contains a putative 
receptor binding domain (amino acid positions 339 to 549) 
with 90 to 96.5% amino acid identities to sequences of ATCC 
and Paris strains. 

Using all seven available HCoV-OC43 genome sequences 
for analysis, the KJK S ratios for the various coding regions 
were calculated (Table 5). The highest KJK S ratios in HCoV- 
OC43 genomes were observed at nsp5 (0.599), followed by HE 
(0.542), S (0.402), and N (0.369), suggesting that these regions 
in HCoV-OC43 are under higher selection pressure. However, 
when using the 29 HK strains (from 2004 to 2011) for analysis, 
the KJK S ratio for the S gene dropped to 0.270. This is likely 
due to the very high KJK S ratio observed at S genes among the 
three genotype A (5.134) strains compared to ratios for geno¬ 
type B (0.332) and genotype C (0.316) strains. 

Estimation of divergence dates. Using the constant popula¬ 
tion size under a relaxed-clock model with an uncorrelated 
exponential distribution, the mean evolutionary rate of HCoV- 
OC43 was estimated at 6.7 X 1CT 4 and 3.6 X 1CT 4 nucleotide 
substitutions per site per year for S and N genes, respectively 
(Fig. 5). These estimates were comparable to previous findings 
of 6.1 X 1CU 4 and 3.6 X 1CU 4 nucleotide substitutions per site 
per year for S and N genes of pigeon herpes encephalomyelitis 
virus, BCoV, and HCoV-OC43 (49). Molecular clock analysis 
using the S gene showed that the tMRCA of all HCoV-OC43 
genotypes was estimated at 1,956.8 (HPDs, 1,940.16 to 
1,966.49) 54 years ago; that of genotype A at 1,964.55 (HPDs, 


1,960.74 to 1,966.73) 47 years ago; that of genotypes B and C 
at 1,984.11 (HPDs, 1968.67 to 1996.31) 27 years ago; that of 
genotype B at 1,996.27 (HPDs, 1,989.83 to 2,001.08) 15 years 
ago; and that of genotype C at 2,001.81 (HPDs, 1,998.66 to 
2,003.82) 9 years ago. Molecular clock analysis using the N 
gene showed that the tMRCA of all HCoV-OC43 genotypes 
was estimated at 1,956.65 (HPDs, 1,933.5 to 1,966.99), also 54 
years ago; that of genotype A at 1,965.58 (HPDs, 1,961.92 to 
1,967) 46 years ago; that of genotypes B and C at 1,989.39 
(HPDs, 1,970.22 to 2,001.91) 22 years ago; that of genotype B 
at 1,997.93 (HPDs, 1,988.21 to 2,002.9) 13 years ago; and that 
of genotype C at 1,998.91 (HPDs, 1,990.87 to 2,003.59) 12 
years ago. These estimates using the two different genes were 
consistent and also were in line with a previous study demon¬ 
strating the tMRCA of HCoV-OC43 strains at 1,944 and 1,957 
by S and N gene analysis, respectively (49). 

DISCUSSION 

The present study represents the first report of possible 
natural recombination among HCoV-OC43 strains, which has 
resulted in the emergence of strains of potentially novel geno¬ 
types. Although HCoV-OC43 was first discovered in 1967 (30), 
genomic studies of HCoV-OC43 have been scarce, with the 
first complete genomes coming from a laboratory strain from 
the ATCC and a clinical isolate, designated Paris, reported in 
2004 (39). This was followed by genomic studies on two HCoV- 
OC43 strains detected in 2003 and 2004 in Belgium, showing 
that the Belgium strains were genetically distinct and that 
HCoV-OC43 could have originated from recent zoonotic 
transmission (47, 48). It also was found that the Paris isolate 
may be cross-contaminated with the ATCC strain, which ex¬ 
plains their close genetic relatedness (48). A later study from 
France analyzing the SI genes of seven HCoV-OC43 strains 
also showed high genetic diversity (43). In this study, we 


FIG. 3. Phylogenetic analysis of nspl to nspl6, NS2a, HE, S, NS5a, E, M, and N genes of seven HCoV-OC43 genomes. The trees were 
constructed by the neighbor-joining method using Kimura’s two-parameter correction, and bootstrap values were calculated from 1,000 trees. A 
total of 738, 1,815, 5,697, 1,488, 909, 861, 267, 591, 330, 411, 46, 2,784, 1,809, 1,563,1,125, 900, 837, 1,244, 4,092, 330, 248, 678, and 1,347 nucleotide 
positions in nspl, nsp2, nsp3, nsp4, nsp5, nsp6, nsp7, nsp8, nsp9, nsplO, nspl2, nspl3, nspl4, nspl5, nspl6, NS2a, HE, S, NS5a, E, M, and N, 
respectively, were included in the analysis. The scale bar indicates the estimated number of substitutions per 50 or 100 nucleotides as indicated. 
The corresponding nucleotide sequences of HCoV-HKUl were used as the outgroups. 
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TABLE 5. Estimation of nonsynonymous and synonymous 
substitution rates in the seven genomes of HCoV-OC43 


Gene 

Function 

HCoV-OC43 substitution 
rates 

K a K s KJK, 

nspl 

Unknown 

0.003 

0.01 

0.24 

nsp2 

Unknown 

0.003 

0.012 

0.245 

nsp3 

Papain-like protease 

0.002 

0.014 

0.17 

nsp4 

Unknown 

0.002 

0.011 

0.142 

nsp5 

3C-like protease 

0.002 

0.004 

0.599 

nsp6 

Unknown 

0 

0.015 

0 

nsp7 

Predicted replicase 

0 

0 


nsp8 

Predicted replicase 

0.002 

0.008 

0.277 

nsp9 

RNA synthesis protein 

0 

0.004 

0 

nsplO 

RNA synthesis protein 

0 

0.003 

0 

nspll 

Unknown 

0 

0.06 

0 

nspl 2 

RdRp 

0.001 

0.009 

0.066 

nspl 3 

Helicase 

0 

0.014 

0.015 

nspl4 

Unknown 

0.001 

0.015 

0.065 

nspl 5 

Endoribonuclease-like 

0.003 

0.015 

0.169 

nspl 6 

Putative methyltransferase 

0 

0.016 

0 

ns2a 

Predicted phosphoesterase 

0 

0.015 

0 

HE 

Helicase 

0.011 

0.02 

0.542 

S 

Spike 

0.014 

0.035 

0.402° 

ns5a 

Unknown 

0.003 

0 


E 

Envelope 

0.003 

0.01 

0.299 

M 

Membrane 

0.004 

0.019 

0.223 

N 

Nucleocapsid 

0.005 

0.013 

0.369 


° The KJK S ratio for S genes of the 29 HK strains from 2004 to 2011 was 0.270. 


showed that there were at least three distinct clusters of 
HCoV-OC43 strains upon RdRp, S, and N gene analysis. One 
cluster, clade A, was formed by the ATCC and Paris strains. 
The other two clusters, clade B and clade C, were formed by 
the present HK strains and Belgium strains BE03 and BE04. 
However, 10 unusual strains displayed incongruent phyloge¬ 
netic positions and belonged to clade B upon RdRp gene 
analysis and to clade C upon S and N gene analysis. These 
results suggested the presence of four different genotypes of 
HCoV-OC43, genotype A (comprising the ATCC and Paris 
strains), genotype B (including Belgium strain BE03 and five 
HK strains from 2004), genotype C (including 15 HK strains 
from 2004 to 2006), and genotype D (including the 10 unusual 
strains: Belgium strain BE04 and 9 HK strains, 1 from 2004 
and 8 from 2008 to 2011). Moreover, genotype D is likely a 
recombinant genotype which has arisen from recombination 
between genotype B and C strains at a region between the 
RdRp and S genes within the genome. To investigate the 
suspected recombination event, complete genome sequences 
of two strains, HK04-01 (genotype C) and HK04-02 (genotype 
D), were determined. Both phylogenetic and bootscan analyses 
showed possible recombination events between genotypes B 
and C in the generation of genotype D strains, a situation 
similar to that reported for HCoV-HKUl (52). The analysis of 
more HCoV-OC43 strains from other countries also will reveal 
the relative prevalence of the different genotypes in different 
localities and the presence of additional genotypes arising from 
other recombination events. 

The recombinant genotype D strains may represent an 
emerging HCoV-OC43 genotype associated with human infec¬ 
tions. The present study, the first molecular epidemiology 
study on HCoV-OC43 infections with clinical characteristics 


presented, revealed genetic evolution into different genotypes 
over time. In a previous study from Belgium, three phyloge¬ 
netic clusters were identified based on S gene analysis, the 
ATCC cluster and two clusters containing four 2003 strains 
and three 2004 strains, respectively, suggesting different tem¬ 
poral patterns among different clusters (48). In this study, 29 
HK strains collected during a 7-year period were included to 
better elucidate the genetic evolution of HCoV-OC43 over 
time. None of the contemporary strains belong to genotype A, 
which consisted only of ATCC and Paris strains that likely were 
isolated 44 years ago. Five of the HK strains from 2004, to¬ 
gether with Belgium strain BE03 from 2003, belonged to ge¬ 
notype B. Fifteen HK strains from 2004 to 2006 belonged to 
genotype C. One HK strain from 2004, eight HK strains from 
2008 to 2011, and Belgium strain BE04 belonged to genotype 
D. While only 1 of the 18 HK strains from 2004 belonged to 
genotype D, all 8 HK strains from 2008 to 2011 belonged to 
this recombinant genotype. This suggests that new genotypes 
of HCoV-OC43 have evolved over time, with the most recent 
HCoV-OC43 strains circulating in our population being dom¬ 
inated by genotype D, which likely has arisen from recombi¬ 
nation as early as 2004. Molecular clock analysis using S and N 
gene sequences suggested that the most recent common an¬ 
cestor of all HCoV-OC43 genotypes emerged in the 1950s 
(mean, 1957), while genotype B and C emerged in the 1980s 
(means, 1984 and 1989 by S and N gene analysis, respectively), 
genotype B emerged in the 1990s (means, 1996 and 1998 by S 
and N gene analysis, respectively), and genotype C emerged in 
the late 1990s to early 2000s (means, 1999 and 2002 by S and 
N gene analysis, respectively). Although the tMRCA of the 
recombinant genotype D strains could not be studied by mo¬ 
lecular clock analysis, the detection of a genotype D HK strain 
and the reported Belgium strain BE04 from 2004 suggested 
that this genotype has emerged no later than that year. More¬ 
over, seven of the eight genotype D HK strains from 2008 to 
2011 were associated with pneumonia, especially in the elderly, 
suggesting that this emerging, recombinant genotype is associ¬ 
ated with more severe disease. However, molecular epidemi¬ 
ology studies involving a larger number of strains and from 
different geographical areas are required to better understand 
the molecular evolution of HCoV-OC43 and the relative 
pathogenicity of the different genotypes. Continuous studies 
also are warranted to detect the emergence of new genotypes 
and recombinants of HCoV-OC43 as well as other human 
coronaviruses and to assess their significance and potential in 
causing future epidemics. Nevertheless, it should noted that 
the amplification and sequencing of a single gene may not be 
sufficient to define the genotypes of HCoV-OC43, HCoV- 
HKU1, HCoV-NL63, and probably other coronaviruses (36, 
52). Given that recombination events are not uncommon 
among these human coronaviruses, the amplification and se¬ 
quencing of at least two gene loci, probably one from ORFlab 
(e.g., RdRp or helicase) and one from HE to N (e.g., S or N), 
should be performed to more accurately understand their mo¬ 
lecular epidemiology and reveal novel genotypes due to recom¬ 
bination events. 

Although MHV is, historically, the most well-studied coro- 
navirus for recombination in in vitro studies (10, 28), there is 
increasing evidence for natural recombination in other coro¬ 
naviruses, some of which lead to the generation of new strains 
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FIG. 5. Estimation of the time to the most recent common ancestor for HCoV-OC43 genotypes. The time-scaled phylogeny was summarized 
from all MCMC phylogenies of the S and N gene data set, which were analyzed under the relaxed-clock model with an uncorrelated exponential 
distribution in BEAST (version 1.6.1). A, genotype A; B, genotype B; C, genotype C. Genotype D strains are in boldface. 
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or genotypes. In feline coronavirus (FCoV), FCoV type II 
strains have been found to have originated from a double 
recombination between FCoV type I and canine coronavirus 
(CCoV) in the M gene (14). Novel CCoV type II strains also 
have been suggested recently to have originated from double 
recombination with porcine transmissible gastroenteritis virus 
in the S gene (6). In our previous study on the natural recom¬ 
bination of HCoV-HKUl in the generation of different geno¬ 
types, the sites of major recombination were localized at nspl6 
just upstream of the stop codon of ORFlab, at the end of nsp6, 
and upstream of nsp5 (52). Recently, we also have demon¬ 
strated natural recombination events among different bat coro- 
naviruses, including those which could have accounted for the 
emergence of civet SARSr-CoV (22). It was found that civet 
SARSr-CoV strain SZ3 was a potential recombinant of horse¬ 
shoe bat SARSr-Rh-BatCoV strains Rp3 from Guangxi Prov¬ 
ince and Rfl from Hubei Province, with a recombination 
breakpoint identified at the nspl6/S intergenic region (22). In 
the study on Rousettus bat coronavirus HKU9, which belongs 
to novel Betacoronavims subgroup D, recombination was iden¬ 
tified at nsp3 and at the nspl5/16 junction (21). In this study, 
possible recombination was observed mainly within ORFlab at 
the nsp2/nsp3 junction, nspl2/nspl3 junction, and NS2a/HE 
junction of HCoV-OC43 genomes. ORFlab is the region most 
susceptible to recombination in the coronavirus genomes. Fur¬ 
ther studies are required to better understand the role of 
recombination in various coronaviruses and the common sites 
of recombination in coronavirus genomes. 

The present results also revealed previously undescribed 
features in the HCoV-OC43 genomes. Although the TRS pre¬ 
ceding most ORFs in HCoV-OC43 genomes has been de¬ 
scribed previously, the putative TRS shared between NS5a and 
E suggested for the ATCC strain (39) was absent from the 
genomes of BE03, BE04, HK04-01, and HK04-02 as a result of 
a 12-nucleotide deletion. While we identified a potential alter¬ 
native TRS, a 39-bp sequence homologous to the leader se¬ 
quence also was found upstream of NS5a. Such a homologous 
sequence has been suggested as a compensation mechanism 
for the absence of TRS in HCoV-NL63 and HCoV-HKUl (34, 
35). Further experiments are required to ascertain the mech¬ 
anism for the translation of NS5a and E protein in HCoV- 
OC43. The analysis of the KJK S ratios for different regions of 
the HCoV-OC43 genome also revealed interesting findings. 
The KJK S ratios were relatively high at nsp5, HE, S, and N 
regions of the HCoV-OC43 genome, suggesting that these 
regions were under higher selective pressure. The nsp5 of 
HCoV-OC43 encodes the chymotrypsin-like protease 3C-like 
protease (3CL pro ) that is important for the proteolytic cleavage 
of the large polyprotein encoded by ORFlab. The 3CL pro of 
SARS-CoV also has been found to induce cellular apoptosis 
(27). Both HE and S are major viral membrane glycoproteins 
which may be important for tissue tropism, cell attachment, 
and eliciting neutralizing antibodies (20, 59, 60). The corona¬ 
virus N protein, a highly phosphorylated protein required for 
viral replication, is also an immunogenic protein which elicits 
subgroup-specific antibody responses (12, 21, 51, 58). Interest¬ 
ingly, the KJK S ratios at the S gene had dropped from 5.134 
among genotype A strains isolated in the 1960s to 0.270 among 
the 29 HK strains from 2004 to 2011. This may reflect the rapid 
evolution of spike protein to adapt to a new host soon after 


interspecies transmission, as HCoV-OC43 was thought to have 
originated from zoonotic transmission, sharing a common an¬ 
cestor with bovine coronavirus dating back to 1890 (47). Fur¬ 
ther studies on more ancient and contemporary strains are 
required to better understand the selective pressure at differ¬ 
ent regions of the HCoV-OC43 genome and its significance in 
terms of protein function and evolution. 
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